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We report on nonlocal transport in superconductor hybrid structures, with ferromagnetic as well as 
normal-metal tunnel junctions attached to the superconductor. In the presence of a strong Zeeman 
splitting of the density of states, both charge and spin imbalance is injected into the superconductor. 
While previous experiments demonstrated spin injection from ferromagnetic electrodes, we show 
that spin imbalance is also created for normal-metal injector contacts. Using the combination of 
ferromagnetic and normal-metal detectors allows us to directly discriminate between charge and 
spin injection, and demonstrate a complete separation of charge and spin imbalance. The relaxation 
length of the spin imbalance is of the order of several /im and is found to increase with a magnetic 
field, but is independent of temperature. We further discuss possible relaxation mechanisms for the 
explanation of the spin relaxation length. 



I. INTRODUCTION 

When a spin-polarized current is injected from a fer- 
romagnet into a spin-degenerate metal, it creates a non- 
equilibrium spin accumulation which is described by a 
difference in the occupation probabilities of states for 
spin up and down. Spin injection from ferromagnets 
has been demonstrated both for normal metals^ and 
spin-degenerate superconductors.^ Recently, spin injec- 
tion from ferromagnets into superconductors in the pres- 
ence of a strong Zeeman splitting has been observedj^j^ In 
this case, spin accumulation mainly originates from the 
spin-dependent density of states in the superconductor;^ 
rather than a difference in occupation probabilities. In 
this article, we extend our previous work and demon- 
strate efficient spin injection from a normal metal into 
a superconductor with a large Zeeman splitting. Using 
both ferromagnetic and normal-metal contacts, we also 
show the separation of spin and charge imbalance in a 
single sample. 

II. SAMPLES AND EXPERIMENT 

Samples were fabricated by e-beam lithography and 
shadow evaporation techniques. In a fist step, a super- 
conducting (S) aluminum wire of thickness ^ai ~ 15 nm 
is created. The aluminum wire is oxydized in situ to 
form a thin but pinhole-free tunnel barrier by exposing 
it to 0.5 Pa of pure oxygen for five minutes. After the ox- 
idation, counterelectrodes of ferromagnetic (F) iron and 
nonmagnetic (N) copper are deposited under a second 
and third angle respectively. Care was taken to have no 
overlaps between the N and F electrodes in the proximity 
of the tunnel contacts. 

We investigated samples with a different number of N 
and F electrodes and variations in the contact distances, 
but otherwise similar parameters. An overview of the 
sample parameters is given in Table HI Sample C, which 




FIG. 1. (color online) Scanning electron microscopy image 
of a section of sample A. An iron (F) and two copper (N) 
electrodes form tunnel contacts with a superconducting alu- 
minum (S) wire. Examples of two measurement configura- 
tions for nonlocal measurements using a normal-metal injector 
and either a normal-metal (NISIN) or ferromagnetic (NISIF) 
detector are indicated. 



has only ferromagnetic junctions, is the same as the sam- 
ple labeled FISIF in Ref. (3|, and is included here only for 
comparison with the new results on samples with mixed 
junction types. Figure [1] shows a scanning electron mi- 
croscopy image of a part of sample A as well as a scheme 
of the measurement setup. The results presented in this 
paper stem from sample A unless explicitly stated. 

All measurements were performed in a dilution refrig- 
erator at temperatures down to T = 50 mK with the 
magnetic field in the plane of the contacts parallel to 
the iron leads, as indicated in Fig. [T] A voltage Vjnj 
consisting of a dc bias and a low-frequency ac excita- 
tion was applied to one tunnel contact, called injector, 
and the ac part of the resulting current /inj was mea- 
sured by standard lock-in techniques to obtain the lo- 
cal conductance gioc = rf^inj/^Mnj- Simultaneously, the 
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ac current /dct through a second contact, called detec- 
tor, was measured to determine the nonlocal conductance 
gni = dIdet/dViny For details of the setup see Ref. @. To 
account for the slight variation in the conductance from 
contact to contact, we plot the normalized nonlocal con- 
ductance ^ni — Sni/GinjGdct throughout the paper, were 
Ginj and Gdct are the normal-state conductances of the 
two junctions. 

The nonlocal conductance was measured for different 
contact configurations, where both injector and detector 
could be either normal (N) or ferromagnetic (F). These 
configurations will be labeled by AISIA, where A and A 
denote the injector and detector contacts, respectively. 
Two possible configurations with a normal-metal injec- 
tor and normal-metal (NISIN) or ferromagnetic (NISIF) 
detector are shown as examples in Fig. [1] The distance 
between injector and detector contact is denoted by d. 



III. MODEL 

Before showing the experimental results, we would like 
to briefly describe the model we have used to analyze our 
data. 

In order to describe the local conductance of the injec- 
tor junctions, we use the theory of tunneling in supercon- 
ductors in high magnetic fieldii— The contribution of a 
single spin projection cr = ±1 to the tunnel conductance 
is given by 



9a 



G 



f (1 



aPinj) / n,{E)fdE, 



(1) 



where Pjnj is the spin polarization of the tunnel conduc- 
tance, n(j[E) is the normalized quasiparticle density of 
states in the superconductor for a single spin, and 
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is the derivative of the Fermi function. The density of 
states na-{E) is calculated by standard methods, includ- 
ing the pair-breaking parameter F, the spin-orbit scatter- 
ing strength 6so, and the Zeeman effect.—"^ The injector 
conductance is given by the sum of the two spin contri- 
butions. 
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whereas the differential spin current is proportional to 
their difference 

ffi - 5t « J [Pini ("i + "t) + ("i - %)] fdE. (4) 

As mentioned in the introduction, two terms contribute 
to spin injection. The the first term is responsible for 
spin injection from a ferromagnet into a spin-degenerate 
superconductor j^i^iiS The second term is due to the spin- 
dependence of the density of states. It appears only in the 



presence of a Zeeman splitting, but does not require a fer- 
romagnet for spin injection. In previous experiments on 
high field spin-injection into superconductors using fer- 
romagnetic j unctions f2,iiL both terms contributed to spin 
injection, while in the present experiment we probe ex- 
clusively the second term. 

The injected quasiparticles create both charge and spin 
imbalance in the superconductor. We describe both on 
an equal footing by a straightforward extension of the 
simple models discussed in Refs. i and [nl The den- 
sities of non-equilibrium charge and spin for each spin 
band are denoted by Q* and So-, respectively. Nonequi- 
librium charge and spin relax over time scales tq* and 
Ts , respectively, leading to an exponential decay over the 
two relaxation lengths Aq* = ^JDtq* and A5 = V Dts, 
where D is the diffusion constant of the superconductor. 
For one-dimensional diffusion along the superconducting 
wire in the geometry of our experiment, their injection 
rates are given by 
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where /* accounts for the fractional quasiparticle 
charge and A is the cross section of the superconduct- 
ing wire. The current flowing out of the detector is given 
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Combining injectioUj relaxation and detection in the 

same way as m Ref. lUj, we obtain 
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where the contribution due to charge imbalance is 
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and the contribution due to spin imbalance is 



5ni = —p, -Pdot-;7-r exp(-d/As) 
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Here d is the distance between the contacts, and pn is 
the normal-state resistivity of the superconductor. These 
two equations form the basis of our data analysis. We 
note that the charge-imbalance signal is always positive, 
whereas the spin-imbalance signal can have either sign. 
Also, detection of the spin signal requires a finite Pdct- 
Therefore, a normal-metal detector will measure only the 
charge signal, whereas a ferromagnetic detector measures 
the sum of the charge and spin signals. Comparing the 
signals of ferromagnetic and nonmagnetic junctions we 
can therefore discriminate the two contributions. 
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TABLE I. Overview of sample properties. Junction properties: number of normal metal (N) and ferromagnetic (F) tunnel 
junctions, range of normal-state tunnel conductances Gn/f, and range of contact distances d. Aluminum film properties: film 
thickness tAi, resistivity pAi, critical temperature Tc, critical magnetic field Be, coherence length ^, diffusion constant D. 
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FIG. 2. (color online) (a) Differential conductance gioc of an 
NIS junction as a function of the bias voltage Vinj for different 
magnetic fields B. (b) Data for i? = 1.0 T together with a 
fit to Eq. ([Sjl and the contributions of the individual spin 
orientations, (c) Calculated spin injection g_i — g^. 



IV. RESULTS 

Figure ^a.) shows the differential conductance of an 
NIS tunnel junction for different applied magnetic fields 
B. At zero field, gioc is negligible for injector bias in the 
subgap region, \V\ < 200 ^V, demonstrating the high 
quality of the pinhole-free tunnel barrier. Above the gap, 
sharp singularities are seen before the conductance drops 
back to its normal state value at high bias. In an ap- 
plied magnetic field, the density of states broadens due 
to pair-breaking, and the Zeeman splitting is clearly seen 
for large fields. ^■^ The data for B = 1.0 T are given in 
Fig. ^h) together with a fit to Eq. Details of the 

fitting procedure are given in Refs. and [3 From the 
fit, the individual contributions ^-t- and to the tun- 
nel conductance are extracted and indicated in Fig. [l{b) 
as dotted and dashed curves, respectively. The differ- 
ence gi — is plotted in Fig. [H^c) . In the bias window 
of the Zeeman splitting around Vinj ~ ±180 fiY, a sin- 
gle spin band dominates conductance, and consequently 
the spin injection has a maximum, g^ — g^ is the only 
bias-dependent quantity which enters Eq. thus the 
nonlocal signal due to spin imbalance is expected to have 
the same shape as depicted here. 

Figure[3ja) and (b) shows the nonlocal conductance for 
the combination of a ferromagnetic and a normal-metal 
contact (the two right-most contacts in Fig. [T|). In panel 
(a), the normal- metal contact was used as injector, and 
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FIG. 3. (color online) (a) Normalized nonlocal differential 
conductance g^i for the closest pair of electrodes with N as 
injector, F as detector (NISIF) and for several magnetic fields, 
(b) ^ni for the same contacts and magnetic fields, but now 
with F as injector and N as detector (FISIN). (c) Difference 
of the signals from (a) and (b) at B = 1.0 T. (d) Comparison 
of NISIF and NISIN data at B = 0.5 T for large contact 
separation d. The lines in (c) and (d) are model predictions 
(see discussion). 



the ferromagnetic contact was used as detector (NISIF), 
whereas in panel (b) their roles are reversed (FISIN). For 
-8 = 0, the nonlocal signal shown in Fig. ^a) is zero for 
bias values below the energy gap of the superconductor 
and increases steeply as soon as the injector bias exceeds 
the energy gap. In an applied magnetic field, additional 
asymmetric peaks appear at voltages near the gap, which 
first grow in height until B « 0.5 T. For higher fields, 
they start to decline, broaden and move inwards. If one 
swaps the role of injector and detector. Fig. ^h), only 
the symmetric part of the signal is found, which quickly 
decreases with applied magnetic field. This signal is due 
to charge imbalance, as has been shown in previous ex- 
periments on NISIN structures j^iii Using Eq. ([7]) we can 
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extract the spin signal alone by subtracting the charge 
signal seen in Fig. from the charge+spin signal in 
Fig.EJa). This is highlighted by Fig. (He) for the data at 
B = 1.0 T. The difference signal shows the asymmetric 
peak structure as well as slight side extrema of different 
sign than the peaks. 

A direct comparison of the charge and spin signals for 
large contact separations d w 6 /im is given in Fig.[3][d). 
Here, the two right-most contacts in Fig. [T]were used as 
detectors, and the same normal-metal contact was used 
for injection. Therefore, the nonequilibrium quasiparticle 
populations probed by the two detectors are essentially 
the same. For these contact distances d, the charge im- 
balance has relaxed and the NISIN signal has thus van- 
ished. In contrast, the asymmetric peak structure of the 
NISIF signal is still visible. This directly probes the spa- 
tial separation of charge and spin imbalance. 

The sohd curves in Fig. ^c) and (d) are calculations 
of the nonlocal signal according to Eq. © and will be 
discussed below. 
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FIG. 4. (color online) (a) Nonlocal conductance g^i for sev- 
eral contact distances d at fixed magnetic field B = 1.0 T. 
(b) Semilogarithmic plot of the peak area A as a, function of 
contact distance d for different magnetic fields B. The lines 
are fits to an exponential decay, (c) Nonlocal conductance gni 
for different temperatures at fixed magnetic field and contact 
distance, (d) Normalized peak area A(T)/A(T = 50 mK) 
for small contact distances d as a function of temperature T, 
shown for all three samples. The line is a guide to the eye. 

In Fig. m^a) the nonlocal conductance ^ni is shown for 
different contact distances d at a fixed magnetic field 
B = 1.0 T. The symmetric high-bias signal due to charge 
imbalance is visible only for the shortest injector-detector 
distance of ci = 0.5 /xm whereas the asymmetric peaks are 
present at the same injector bias value for all distances 



and are visible even for the largest distance of d = 6.4 fim. 

To further analyze the dependence of the signal on con- 
tact distance, we extracted the charge signal by calcu- 
lating the symmetric part gs = [5ni(T4nj) + 5ni(-^nj)] /2, 
and the spin signal by calculating the antisymmetric part 
5a = [5ni(V^nj) ~ gni(-Vinj)] /2. For the spin signal, we 
then calculated the peak area A by integrating nu- 
merically from Vinj = —290fiV to T^nj = O^V. 

Figure HJ^b) shows the peak area A as a function of 
contact distance d for different magnetic fields B on a 
semilogarithmic scale. For better visibility, the datasets 
have been offset vertically. The data can be fit to an 
exponential decay, shown by the lines in Fig. Iljb), ex- 
cept for the closest contact distance, which was there- 
fore excluded from the fit. From these fits the spin- 
imbalance relaxation length is extracted. To extract 
the charge-imbalance relaxation length Aq» , we also fit 
(js Sit Vinj = 520 fiY to an exponential decay (not shown). 
The results of these fits are shown in Fig. [S] 

In Fig. \^c) we show the nonlocal conductance ^ni for a 
fixed magnetic field of i? = 1.0 T and different tempera- 
tures T. The infiuence of temperature on the spin signal 
is to decrease the peak height and broaden the peaks. 
However, the influence is not just a thermal broadening, 
as the peak area decreases significantly when the tem- 
perature is increased. To analyze the decrease in greater 
detail, we normalized the peak area A to the value at 
T = 50 mK and found a similar decrease with tempera- 
ture for all three samples as can be seen in Fig. Sfd). 

The relaxation lengths As and Aq. obtained from the 
exponential fits described above are shown in Fig. [SJa) 
as a function of magnetic field B for all three samples. 
Xs increases monotonously with B, but one observes for 
all three samples a change in slope at B « 0.5 T, which 
is most pronounced for sample B. For the highest fields 
close to Be, the behavior gets less systematic. Since the 
absolute values of the peak areas decrease significantly 
when approaching Be, the reduced signal-to- noise ratio 
serves as a plausible explanation. Ag., in contrast, is 
largest for B = and reduces quickly in a magnetic 
field. While As show a pronounced dependence on the 
magnetic field B, no significant change with temperature 
T can be found, as it is shown for all three samples in 
Fig.ETb). 

Figures [Sfc) and (d) present the charge-imbalance and 
spin-imbalance relaxation times tq. and ts- These were 
determined from the relaxation lengths shown in Fig.[5l^a) 
using the relations rg. = Ag, /D and rs = A| /D and the 
known diffusion coefficients D given in Table HI 

At zero magnetic field, the charge imbalance relaxation 
rate Tq} is determined by the combined effect of inelas- 
tic scattering and elastic impurity scattering in conjunc- 
tion with the gap anisotropy of aluminum. ^§ In addi- 
tion, magnetic pair-breaking perturbations increase the 
relaxation rate.iS Previous experiments on charge imbal- 
ance have shown that these two contributions are addi- 
tive at low temperature, i.e., Tq} = Tq"^ + where Tq"^ 
denotes the relaxation rate at zero field, F is the pair 



5 



- {b) O sample A, B= 1.00 T 

□ sample B, B= 0.50 T 
7.0 - O sample C, B- 0.75 T 




2HgB / kT 
10 20 30 40 50 



(^) ♦ sample A 
• • sample B 

■ sample C 



O sample A 
O sample B 
□ sample C 



0.2 0.4 0.5 O.S 1.0 

B (T) 



0.0 0.5 1.0 1.5 

B (T) 



FIG. 5. (color online) (a) Spin (open symbols) and charge 
imbalance (closed symbols) relaxation lengths Xq* and As as 
a function of the magnetic field B. (b) Spin relaxation length 
As as a function of temperature T at fixed magnetic field B. 
(c) Charge imbalance relaxation time rg* as a function of 
the magnetic field B, together with fits (see text), (d) Spin 
relaxation time rs as a function of the magnetic field B, with 
an attempt to fit with exp(2/iB-B/fcBT). 



breaking parameter, and /3 is a coherence factor iii Us- 
ing r = [B / /2 for a thin film in parahel magnetic 
field, we obtain tq* = (tq~^ + hB^)~^. Fits to this de- 
pendence are shown in Fig. [SJc). As can be seen, tq* is 
a few nanoseconds at zero field, and then rapidly drops 
at higher fields. 



Figures [Sfd) shows the spin relaxation time ts as a 
function of magnetic field i?, where we also show the 
ratio of Zeeman splitting and temperature, 2^bB /k^T , 
at the top of the panel, ts is a few nanoseconds at 
small fields, similar to tq, . At higher fields, it in- 
creases to 10 — 20 ns, and exceeds tq* by at least two 
orders of magnitude. When the Zeeman splitting is 
not much larger than the thermal broadening, we ex- 
pect S^/S^ oc exp(2//B^/^Br). Therefore, the relax- 
ation time might show a similar dependence on magnetic 
field, as assumed in Ref. d. To check this assumption, 
we have plotted an attempt to fit ts cc eyip{2^^B / k^T) 
in Fig. md). As can be seen, the fit is possible at best 
in small fields, but clearly fails to describe the high-field 
data. 



V. DISCUSSION 

In our previous work^i^ we reported on spin injection 
from ferromagnetic contacts into a superconductor in the 
presence of a Zeeman splitting. In that case, both the 
finite spin polarization of the injector, Pmj, and the spin- 
dependence of the density of states, nj^ — n-|-, contribute 
to spin injection, as inferred from Eq. — n-f- is 

responsible for the asymmetric shape of the signals. As 
explained in Ref. the role of Pinj is only to increase 
the amplitude of the positive peak, and to decrease the 
amplitude of the negative peak (in other words, there is 
an overall upward shift of the signal). As can be seen 
in Figs. [3] and m the amplitudes of the peaks at positive 
and negative bias are the same in the present experiment, 
where spin is injected from a normal metal with Pjnj = 0. 

For a quantitative comparison with the model, the non- 
local signals calculated from Eq. ^ for a spin-degenerate 
injector are plotted as solid lines in Figs, and (d). 
The factors {g^ — g^)/Gi-a\ and Pdet are obtained from 
the fits of the local conductance of the injector and de- 
tector junctions, as shown for example in Fig.[5Jb). The 
relaxation length As is obtained from the exponential 
fits in Fig. [D^b). The normal-state resistance per length 
of the aluminum wire is known from measurements at 
T = 4.2 K, so that we can calculate the factor pj<iXs/2A. 
Thus the nonlocal signal g^^ predicted by our simple 
model can be calculated without free fitting parameters. 
The shapes of the calculated and measured signals agree 
qualitatively, whereas the amplitude of the calculated sig- 
nal is too small. Since the spin injection and detection 
factors as well as the normal-state properties are known 
quite accurately, we suspect that the assumption of a 
single exponential decay length A5 independent of en- 
ergy and bias conditions is the most likely culprit for the 
disagreement. To elucidate this further, we would like to 
discuss some possible relaxation mechanisms of the spin 
signal. 

Nonequilibrium quasiparticles in superconductors 
are subject to several different scattering processes. 
Electron-phonon scattering leads to energy relaxation, 
charge-imbalance relaxation, as well as recombination 
of quasiparticles to Cooper pairsJ^ Spin-orbit scatter- 
ing leads to elastic spin flips, and is therefore expected 
to relax spin imbalanccj^^ It also modifies the density 
of states in the presence of Zeeman splitting.'' Magnetic 
pair-breaking perturbations lead to additional charge- 
imbalance relaxation, and may affect all other scatter- 
ing mechanism by changing the density of states and co- 
herence factors. Magnetic impurities in particular also 
lead to spin flips The spin accumulation in a sin- 
gle spin band, So-, is proportional to the product of 
the density of states and the occupation probabil- 
ity fa-. Therefore, for the net spin accumulation we have 
Si- Sf(x nifi - n^f^. 

Without Zeeman splitting, we have nj, = n-j-, and 
spin accumulation is due to the difference in occupation 
probability — /-[-. Several experiments have probed 
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FIG. 6. (color online) Schematic view of the spin-resolved 
density of states, including Zeeman splitting and spin- 
orbit scattering, and various possible relajcation mechanisms. 
States are partially occupied up to i5 = —eVinj (see discus- 
sion for details), (a) Spin-flip scattering in a spin-degenerate 
superconductor, (b) Proposed two-stage relaxation mecha- 
nism involving spin flips and recombination for bias voltages 
in the range of the Zeeman splitting, (c) Inelastic spin flips to 
the upper Zeeman band, (d) Two-stage spin-flip and energy 
relaxation process at high bias. 



spin injection from ferromagnets into spin-degenerate 
superconductors. In this regime, spin relaxation is 
due elastic spin flips by spin-orbit scattering or magnetic- 
impurity scatteringfi^ as sketched in Fig. ^a). The 
spin-relaxation length in the superconducting state is 
expected to be either the same (spin-orbit scattering) 
or shorter (magnetic-impurity scattering) than in the 
normal-state. We have determined the spin-relaxation 
length in the normal state, An, by nonlocal spin- valve 
experiments for sample C, and found An = 370 nm4 For 
samples A and B, we could not determine An, since these 
samples have only one or two ferromagnetic junctions. 
Hovifever, the aluminum film parameters are similar, and 
we can assume An ^ 400 nm for all three samples. 

For spin injection from a normal metal in the presence 
of a Zeeman splitting, we expect /j, = /t-, and spin ac- 
cumulation is solely caused by the difference in density 
of states, nj. — n-|-. By detailed balance, the net spin- flip 
rate due to either spin-orbit scattering or magnetic im- 
purity scattering should therefore be zero, and not lead 
to any spin relaxation at all. Neither energy nor charge 
relaxation have any effect on spin accumulation. Recom- 
bination reduces the overall number of quasiparticles, but 
does not lead to spin relaxation, as it removes one quasi- 
particle from each spin band. In the energy window of 
the Zeeman splitting, however, recombination will de- 
plete the spin-up band much faster than the spin-down 



band due to the different density of states. Therefore, it 
will indirectly enable spin flips. One possible spin relax- 
ation mechanism in our experiment is therefore a two- 
stage process of recombination and spin-flip scattering, 
as shown schematically in Fig. IH^b). Since the normal- 
state spin-diffusion length An is much shorter than the 
observed As, we can assume that recombination is the 
bottleneck for this mechanism. The impact of spin-orbit 
scattering on recombination in high fields was considered 
theoretically for SIS tunnel junctionsJ^ 

The two stage relaxation mechanism can explain why 
As is very large, since recombination is expected to be 
very slow at low temperatures. However, once (IbB ^ 
k-^T, a further increase of the field will essentially no 
longer change the density of states in the energy range 
E < — eVinj occupied by quasiparticles. Therefore, it 
is not clear why As should continue to increase at high 
magnetic fields. One possible explanation would be in- 
elastic spin flips^" to the upper Zeeman band, as shown 
in Fig. [Ifc). The inelastic spin- flip rate should de- 
crease upon increasing field, since larger energy transfer 
is needed. As second possible explanation is the phonon 
bottleneck of recombination)^ For higher fields, the over- 
all number of nonequilibrium quasiparticles in the energy 
window of the Zeeman splitting increases, and within the 
two-stage relaxation mechanism sketched above, more re- 
combination phonons would be created. These would in 
turn reduce the net recombination rate. 

Also, as can be seen in Fig. |3Jd), the measured spin 
signal extends to higher bias than expected from the spin- 
injection factor — obtained from fitting the injector 
conductance spectra. For bias beyond the energy win- 
dow of the Zeeman splitting, an almost equal number of 
quasiparticles with both spins is injected. Cooling due 
to inelastic scattering will then be more efficient for the 
spin-down quasiparticles, since a larger density of states 
is available at low energy. This will free spin-down states 
at high energy, and might therefore lead to a net spin-flip 
scattering rate from the spin-up to the spin-down band. 
This would have the counterintuitive effect that spin-flip 
scattering can lead to an increase of spin accumulation, 
as sketched in Fig. [SJ^d) , and therefore explain that the 
spin signal extends to higher bias than expected. 

Finally, we would like to address the temperature de- 
pendence. As shown in Fig. [Sfd), the relaxation length 
As is almost independent of temperature from 50 mK to 
500 mK. In our experiment, we have both eVinj ^ A:bT 
and A 3> k^T. Therefore, the nonequilibrium quasi- 
particle distribution is determined mostly by bias, and 
is almost independent of temperature. In addition A 
is almost constant in this temperature interval, so that 
also the density of states and coherence factors are 
almost constant. This may explain the temperature- 
independence of As. A similar behavior was found for 
Aq- . - ^ In contrast to As, the peak area of the spin signal 
shown in Fig.HJd) decreases by about 30 % upon increas- 
ing the temperature from 50 mK to 500 mK. This can 
not be accounted for by a change in the injection fac- 
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tor gj. — g-f, since thermal broadening will not lead to a 
change in peak area. 

It is obvious from this discussion that spin relaxation 
in superconductors in high magnetic fields is a complex 
process, and a detailed quantitative model is beyond the 
scope of this article. 



VI. CONCLUSION 

We have shown spin injection and transport in meso- 
scopic superconductors in the regime of large Zeeman 
splitting, and investigated in detail the role of spin- 
polarized and spin-degenerate injector and detector junc- 



tions. We have found that spin injection is possible from 
a normal metal, whereas a ferromagnet is needed as de- 
tector to observe spin accumulation. For spin injection 
from a normal metal, spin accumulation is purely due to 
the spin-dependent density of states in the superconduc- 
tor in high magnetic fields. Comparing the nonlocal con- 
ductance probed by spin-degenerate and spin-polarized 
detectors, we were able to directly discriminate charge 
and spin imbalance. The spin relaxation length increases 
strongly in a magnetic field, but is found to be almost 
independent of temperature. A detailed explanation of 
the relaxation mechanisms remains an open question to 
theory. 

This work was partially supported by the Deutsche 
Forschungsgemeinschaft under grant BE-4422/1-1. 
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